Characterization of nanostructured material images using fractal descriptors 
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This work presents a methodology to the morphology analysis and characterization of nanostruc- 
tured material images acquired from FEG-SEM (Field Emission Gun-Scanning Electron Microscopy) 
technique. The metrics were extracted from the image texture (mathematical surface) by the vol- 
umetric fractal descriptors, a methodology based on the Bouligand-Minkowski fractal dimension, 
which considers the properties of the Minkowski dilation of the surface points. An experiment with 
galvanostatic anodic titanium oxide samples prepared in oxalyc acid solution using different condi- 
tions of applied current, oxalyc acid concentration and solution temperature was performed. The 
results demonstrate that the approach is capable of characterizing complex morphology character- 
istics such as those present in the anodic titanium oxide. 



I. INTRODUCTION 

The morphology analysis of solid samples is an impor- 
tant research area in Materials Science to characterize 
some of their properties |13] |16l [171 EI] • Generally, the 
technique employed in this kind of application consists 
in the quantitative analysis of the micrographs obtained 
using one (or more) of following techniques, such as 
FEG-SEM (Field Emission Gun-Scanning Electron Mi- 
croscopy), AFM (Atomic Force Microscopy), STM (Scan- 
ning Tunneling Microscopy) or TEM (Transmission Elec- 
tron Microscopy). In any case, the results are a matrix of 
values which expresses the topography of the measured 
sample. 

When the morphology of the material is investigated, 
we observe that each image presents a specific distribu- 
tion pattern. This distribution is quite similar to that 
found in textures classically studied with image anal- 
ysis tools, e.g., the so-called texture analysis methods 
[12]. Therefore, the application of this procedure over, 
for example, FEG-SEM images, could make possible the 
modeling of physical important properties of the mate- 
rials, such as roughness, grain boundaries, morphologic 
defects and total surface area. The quantitative analy- 
sis of nanostructured materials for well behaved samples 
has been described using simple functions of existing soft- 
ware packages. One example is the automatic counting of 
pores in self-organized anodic porous alumina [16 using 
ImageJ [Ij and Gwyddion [18]. The problem is different 
when studying complex morphology characteristics of the 
samples as those ones where the distinction among the 
patterns is not obvious. In this case, a more sophist i- 
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cated analysis must be used and these ones generally are 
not included in those software packages described above. 

Among the methods described in the literature, those 
ones based on fractal analysis have presented excel- 
lent performance in the investigation of complex tex- 
tures, mainly on those synthesized and natural samples 
[H [3l [71 [19] . Actually, nature is rich in self-similar pat- 
terns, that is, structures which repeat themselves under 
different scales. From a mathematical point of view, this 
is also an intrinsic property of fractal objects. There- 
fore, fractal geometry is appropriate to measure such 
kind of structures and, as consequence, self-similarity also 
measures complexity (meaning the level of details along 
scales) which is directly related to spatial occupation in 
the structure. 

Notwithstand that fractal dimension provides a good 
solution in many object identification problems, it is lim- 
ited in the representation of some classes of objects [8] 
due to two main aspects: i) fractal dimension is a real 
number which is insufficient to characterize such objects, 
ii) there are samples that have different patterns but 
present the same fractal dimension. In fact, these ob- 
jects can have a more or less self-similar aspect depend- 
ing on the scale observed [14]. In the literature, differ- 
ent propositions have been presented to solve the above 
drawbacks. The main ones are multifractal [11 , multi- 
scale fractal dimension [15] and fractal descriptors 0. 
Considering that several papers have demonstrated the 
superior perfomance of fractal descriptors dealing with 
texture images over the other ones, we are focused here 
on such approach as a tool for texture discrimination 

We apply Volumetric Minkowski descriptors method- 
ology, which was initially developed in [2 . It is derived 
from the Bouligand-Minkowski fractal dimension. The 
descriptors are obtained by mapping the original gray 
level image (of FEG-SEM data, in this case) onto a three- 
dimensional mathematical surface. Thus, such surface is 
dilated by Bouligand-Minkowski method using spheres 
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with predefined radii. The fractal descriptors are then 
estimated from the volume of dilation for each sphere ra- 
dius. With the growing of the dilation radius, the spheres 
start to interfere among themselves, forming a wavefront 
which is tightly related to the structure of the material. 
It is important to stress out that the dilation process 
captures the arrangement of the topography [2 . Thus, 
these descriptors are capable of providing very rich infor- 
mation about the morphology of the material and, conse- 
quently, are a strong method for nanostructured material 
characterization task. The use of the Volumetric Fractal 
Descriptors applied to nano structured surfaces was ini- 
tially proposed in [9 . In that seminal project, was sug- 
gested that Fractal Volumetric Descriptors could be used 
to characterize and analyze such nanostructures, show- 
ing the discrimination power on two distinct conditions. 
In the present work, an experiment with galvanostatic 
anodic titanium oxide samples prepared in oxalyc acid 
solution using eight different conditions of applied cur- 
rent, oxalyc acid concentration and solution temperature 
was performed, and demonstrate that the proposed tech- 
nique is capable to identify the nanosurfaces. The nature 
of the materials surfaces and its images are a difficult 
problem of image analysis, and the proposed technique 
demonstrate be suitable to characterize and to identify 
nanostructured surfaces. 



II. MATERIALS AND METHODS 

A. Materials 

The samples used in this work were galvanostatic an- 
odic titanium oxide ones prepared in oxalyc acid solution. 
In this electrochemical preparation method, a titanium 
plate is the anode in a two electrode electrochemical cell. 
A platinum plate was used as cathode. Then, the an- 
ode is polarized under constant current condition and an 
oxide film starts to form over the anode following the 
equation: Ti + 2H2O Ti02 + 4i^+ + 4e" It is impor- 
tant to stress that Ti02 is formed by the direct reaction 
between the metal and water over the metal. The surface 
morphology of the oxide is sensible to the experimental 
conditions used. In the present case, it were used differ- 
ent values of applied current, oxalyc acid concentration 
and solution temperature as described in Table [l| Using 
a 2^ factorial design, it were performed 8 titanium oxide 
anodizations, generating, therefore, 8 classes of samples. 
From each class, it were acquired images from 8 differ- 
ent regions on their surface. Therefore, we have a total 
of 64 samples to be used in the model building. Each 
sample is a rectangular piece of the plate, which is mea- 
sured through SEM-FEG technique generating a matrix 
(image) with a resolution of 3072x2060 pixels. Figure [l] 
shows one image per class, illustrating the general aspect 
of the dataset. 



III. RESULTS 

A. Fractal Theory 

Fractals are objects formally defined as a set of points 
whose Hausdorff-Besicovitch dimension (see the concept 
below) exceeds strictly the Euclidean dimension. In prac- 
tice, it is an object generated through a dynamic system 
that presents infinite complexity and self-similarity [14]. 
Here, complexity states for the level of details under dif- 
ferent scales. Self-similarity expresses the fact that if 
one takes the fractal under different observation scales, 
he will observe the repetition of patterns only changed 
by simple geometrical transforms (affine transforms). It 
is important to stress out that these properties are ob- 
served at infinite scale. In nature, we can find a lot of 
objects which present self-similarity and complexity at 
certain levels interval. Then, this is a strong motivation 
to approximate these structures through fractal metric. 
The most used of such fractal metrics is the one called 
fractal dimension. 



1. Fractal Dimension 

The first concept presented in the literature of fractal 
dimension [14 is the same one as Hausdorff dimension, 
which is based on the Hausdorff measure S^^ concept. Let 
A be a subset of and s and S non-negative real values. 
Then: 

i3|(A) = inf{^ ||/7i||"such that {[/J is a d-cover of A}, 
where {Ui} is a J-cover of A if A C Ui^i being < 

\m<s. 

The Hausdorff s-measure io* is given by: 

An interesting and important characteristic of this 
measure is that S^^{A) is always for any s < dn and 
00 for any s > dn- The real value dn is the so-called 
Hausdorff dimension of A, that is: 

dH{X) = inf {s\H'{X) = 0} = sup{s|i^"(X) = 00} . 

A special case of fractal dimension definition is the 
Bouligand-Minkowski (BM) dimension, decribed in the 
following section. 

2. Bouligand-Minkowski 

As well as in the Hausdorff dimension, BM also has 
an associated measure, which is, in this case, an upper 
measure ^ and a lower one Qr defined through: 

^(X, R) = lim inf (X, r) , 
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Figure 1. Titanium oxide dataset. Each image number corresponds to the same number of the experiments column presented 
in Table H 
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Table L Experimental Matrix for the anodic titanium oxide samples preparation using a factorial design procedure 
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(X, R) = lim sup Qr (X, r) , 



where 



qr{X, R,r) 



V{dX © rR) 



with — oo < r < oo, r < 0, 9X is the boundary of X 
and denotes the morphological dilation by an element 
R with radius r. 

The upper and lower dimension, Db and Db respec- 
tively, are defined by: 

S^(X,i?) = inf{r|^(X,i?) = 0}, 



represented in FEG samples using the above Bouligand- 
Minkowski method due to its precision O |l9l [20 . We 
applied a neighborhood approach once the FEG image is 
described in a discrete space. 

The most intuitive way of calculating fractal dimen- 
sion of a FEG image is to map the data onto a 3D gray 
intensity surface. The Figure |2] shows an example of 
the image presented as a 3D gray surface. This is per- 
formed in a simple manner by representing the image 
/ G [1 : M] X [1 : TV] ^ 3? in the surface 5' 



S = {i,jJ{i,j)\{i,j)e[l:M]x[l:N]}, 



(1) 



and 



Db{X, R) = inf{r|2^(X, R) = 0}. 

In a discrete space, like that of digital images here an- 
alyzed, the direct application of above equations is not 
viable. In such situations, a common practice is to em- 
ploy neighborhood techniques. Particularly, here, we are 
interested in the estimation of dimension in 3?^. In this 
space, the BM dimension may be presented through: 

I,.(-Y)=li„.(3-«2^y^), 

^ ^ r^o^ logr ^ 

where V is the volume of the dilated structure and Yr is 
an Euclidean sphere with radius r. 



3. Fractal Descriptors 

As described in the introduction section, the fractal 
dimension is a scalar value and it is not enough to char- 
acterize such complex structures as those presented in 
Figure [l] Therefore, we have developed an improved 
concept related to fractal dimension which is not a real 
number but a vector [8 . This mathematical object was 
called fractal descriptor and its fundamental aspects is 
described in the following paragraphs [3l [3 [121 US • 

We propose to extract the morphological properties 
from the samples through an analysis based on fractal 
geometry of the FEG images. In this way, an initial pro- 
cedure is to estimate the fractal dimension of the data 



f{ij) = {l,2,...,max_m/iie}|/ = (2) 
being max.value the maximum value in the FEG data. 
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Figure 2. Texture image mapped onto a 3D surface, (a) Plain 
image, (b) Surface, obtained considering the gray levels as z 
axis. 

A classical way of calculating the BM dimension of a 
surface in discrete space is dilating it with spheres varying 
the radius r, which is presented in Figure [3) This Figure 
illustrates the dilation process for two different values of 
r. Using this approach, we compute the dilation volume 
V{r) (number of points inside the dilated surface) for 
each value of r. The dimension is given through: 



= 3 — lim 



log(r) 
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For computational efficiency purposes, an usual solution 
to calculate V{r) is to employ the exact Euclidean Dis- 
tance Transform (EDT) [6], which is for a surface S in 
three dimensions given by: 

EDT{p) = mm{d{p, q)\q & U\S}, (3) 

where d is the Euclidean distance and U\S corresponds 
to the complement of S taken over a cube U (universal 
set) which contains S. When we are dealing with exact 
EDT, the distance has predefined values: 

^ = 0,1,V^,..., (4) 

where 

l(=D = {d\d={f+fy^^;i,j€N} (5) 

Initially, we define the set gr{S) of points at a distance 
r from S: 

gr{S) = {{x,y,zmx - + {y - Syf + {z - S,f]'/^ 

(6) 

where Sx^ Sy^ Sz are the coordinates of points in S. 
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Figure 3. Dilated surfaces with different radii, (a) Original 
surface, (b) Radius 2. (c) Radius 10. 



The dilation volume V{r) is given by: 

nr) = ^Q(i), (7) 

i=l 

where Q{i) is provided through the following expression: 

{x,y,z)eU 

where x states for the indicator function. 

Therefore, V{r) and r present a power law relation. 
The descriptors themselves are obtained using the fol- 
lowing relation u: 

u : log(r) log(V(r)), 

where V acts as a fractality measure and r is the scale 
parameter. 

IV. DISCUSSION 

Several times, to extract the information from function 
ii, it is necessary a mathematical space transformation, 
such as Fourier or Principal Component transform. 

Here, we apply the Principal Component (PC) trans- 
form [5 of the fractal descriptors in the classification of 
the different investigated materials. Figure [4] illustrates 
the ability of discrimination through BM fractal descrip- 
tors. 

The performance of the proposed approach is verified 
in a task of classification of titanium oxide films prepared 
using different experimental conditions as described in 
the Section [lI A| Each condition corresponds to one class. 

The fractal descriptors are extracted from the images 
and the Principal Component Analysis (PCA) is applied 
over each set of fractal descriptors. The components are 
thus classified by Linear Discriminant Analysis (LDA) 
0, using leave-one-out cross-validation procedure. To 
check the quality of the results, they were compared to 
two classical texture analysis methods, that is, Fourier 
[To] and histogram entropy [10 . To guarantee that the 
comparison is valid, the same PCA-LDA methodology 
was applied also with these last methods. 

Figure [5] illustrates the behavior of PC components rel- 
ative to their correctness rate in the classification process 
for each methods used. It shows the number of compo- 
nents up to 16 to detect the necessary number of principal 
components to be used. 

We observe that the rate of correctness increases as the 
number of components increases and the Fractal descrip- 
tors achieve its maximum values faster than the other 
methods used. In this case, it is necessary to use 14 
components. The results presented in the PC graph are 
summarized in the Table |n| showing the best result ob- 
tained by each descriptor, with the measure error and the 
number of components used to achieve such result. After 
the maximum correctness, the performance decrease for 
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Figure 4. Discrimination performance of proposed descriptors. 



every method which means that the use of more compo- 
nents in the LDA procedure only add noise to the set of 
features, damaging the performance of the descriptors. 

Now, we see that fractal descriptors showed the best 
correctness rate result with a 8% of advantage over the 
second best approach, Fourier descriptors. This result 
was expected due to the intrinsic ability of fractal de- 
scriptors in extracting a rich information of the topogra- 



Method Correctness Rate (%) Number of descriptors 



Entropy 
Fourier 
Fractal 



59.38 =b 0.01 
81.25 =b 0.01 
87.50 ± 0.01 



11 
15 
14 



Table II. Correctness rate for each compared descriptor. 
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Figure 5. Correctness rate according to the number of PC 
components in each compared method. 



V. CONCLUSIONS 

This work showed an apphcation of fractal descriptors 
to the classification of samples of titanium oxide material 
under different experimental conditions. 

We compared the performance of fractal to other clas- 
sical texture descriptors approaches in the literature, 
to know, histogram entropy and Fourier. The results 
showed that fractal approach obtained the best result, 
providing the most accurate classification. We verified, 
as expected, that fractal geometry is a powerful tool to 
describe such nanoscale image. This is explained by the 
flexibility of fractals in modeling topographies arising 
from natural systems, like that present in titanium oxide 
samples. 

The result encourages the research and enhancing of 
novel fractal-based approaches, applied to many challeng- 
ing problems related to the discrimination and descrip- 
tion of materials under different experimental conditions. 

It is important to stress out that this procedure can be 
used also for different natural or synthetic samples which 
can be found in the literature from different areas of the 
knowledge. 



phy of the surface. Such topography is directly related 
to important physical aspects of the material, like rough- 
ness, grain boundaries, morphologic defects, reactivity 
and total surface area. In this sense, we observe that the 
proposed descriptors are capable of accurately capturing 
nuances which are essential in the discrimination skill of 
our vision system. It is important to notice, however, 
that the use of descriptors associated to the classifier 
computational method turns possible the analysis of de- 
tails with a complex processing framework which ensures 
a more precise and robust solution for the discrimination 
problem. 
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